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ABSTRACT 

Three simultaneous rocket measurements of D region 
ionization sources and electron and ion densities have been 
made in one day. Electron density profiles for solar zenith 
angles of 53.2, 27.8 and 48.3° are displayed. The profiles 
ar6 found to be quite similar to those reported for mid- 
latitudes under quiet conditions. Using simultaneous 
measurement of the X-ray flux together with NO ionization 
estimates based on recent measurements of NO by Meira (1971) 
and OoC^A ) ionization rates by Huffman et al (1971) , the 

A g 

ion pair production function is computed. It is demonstrated 
that the charged particle distribution can be explained using 
this production function and an effective recombination 
coefficient compatible with recent laboratory measurements 
of the dissociative recombination coefficient for electrons 
with hydrated cluster ions. 



INTRODUCTION 


The D region is the least understood part of the iono- 
sphere. Although several sources of ionization have been 
proposed, uncertainties of measurement have led to wid? 
speculation concerning their relative importance. Lati- 
tudinal, seasonal, and diurnal behavior of the most easily 
measured parameter, the electron density, are ill -defined. 

The species of ions present are either unknown or at best 
known only qualitatively. It is therefore, necessary to 
improve measuring techniques and to conduct experiments in 
which the sources of ionization, the resultant electron ard 
ion density distribution, and the species of ion are all 
determined simultaneously. 

Three Nike-Apache rockets were launched from the 
Thumba Equatorial Rocket Range in India, near the geomagnetic 
equator, on March 19, 1970. The instrumentation on each 
rocket together with the local time of launch and solar 
zenith angle are listed in Table 1. Additional details 
relating to the flight and the ion composition measurements 
are dealt with in Goldberg and Aikin (1971) . The present 
paper is concerned with the electron and positive ion 
densities and their relationship to the different sources 
of ionization. A production function is derived and used 
in conjunction with the charged particle density to evaluate 
the effective recombination coefficient, or. 

ELECTRON DENSITY EXPERIMENT 

The electron density distribution has been determined 
from a radio wave absorption experiment developed by Kane 
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(1969). The output of a CW transmitter is alternately fed 
to antennas oriented magnetic east-west (ordinary mode) and 
north-south (extraordinary mode) . The amplitudes of the 
resulting signals are monitored by a rocket-borne receiver 
with a linearly polarized antenna. 

The azimuth of the rocket trajectories nearly coincided 
with magnetic west to permit radio wave transmission quasi- 
transverse to the Earth's magnetic field. The frequencies 
3030 and 1865 KHz were chosen for the experiment on each 
of the three payloads. The switching of the propagation 
mode was accomplished every 0.5 seconds or approximately 
once each 2.5 to 4 rolls of the payload. Figures la, lb. 
and lc give the resultant profiles together with their 
respective scatter of data points for the extraordinary 
wave mode on each of the three payloads. Ordinary wave 
data contain more scatter than extraordinary and appear 
to be less accurate. All datavrere reduced using total 
absorption rather than differential absorption techniques 
for increased accuracy in the results. 

ION DENSITY PROFILES 

Total positive ion density was measured in situ usiug 
the Gerdien condenser technique described by Bourdeau et al., 
(1966). A schematic of the instrument used in the present 
experiments is shown in Figure 2. On both flights 14.425 
and 14.424, the Gerdien was operated with the collector (B) 
at a fixed D.C. bias of -4 volts with respect to the grounded 
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payload housing. Assuming total current collection for 
the flow passing through the instrument, the ion density 
profile obtained for 14.424 is illustrated in Figure lb. 

Above 75 km there is excellent agreement between the positive 
ion and electron density. There is an increasing deviation 
between the two profiles below 70 km which is attributed to 
the presence of negative ions. However, the exact concentra- 
tion of these negative ions cannot be determined because 
of inaccuracies which may exist in the radio absorption 
data and the Gerdien flow. 

The positive ion density data for rocket 14.425 was not 
in good quantative agreement with the electron density pro- 
file and it was necessary to normalize the data to the 
electron density profile shown in Figure la by shifting it 
as indicated. The deviation between Gerdien and radio 
absorption data on this flight is believed to be primarily 
due to Gerdien data interpretation. This is caused 
by oblique shock effects which pinch off the effective flow 
cross-section, i.e. effective aperture area, leading to 
erroneously small positive ion densities. Also, NASA 14.425 
traveled at a velocity of Mach 2.2 in the region of discrepancy, 
whereas the better agreement of 14.424 occurred at Mach 3.2. 
These differences and their effects on pinching the flow are 
currently under investigation. 

QUALITATIVE RESULTS 

The data discussed in the previous two sections have 
been combined to form the three profiles for electron 
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density shown in Figure 3. These data were taken under 
quiet time conditions with Ap 10 5. 

The profiles exhibit the shape characteristics of the 
lower ionosphere at midlatitudes. There is a Chapman like 
D region centered between 75 and 85 km. Above 85 km there 
is a steep gradient of electron density which characterizes 
the bottom of the E region. 

The most interesting aspects of the profiles are the 
changes which occur with solar zenith angle (x) . The E 
region is seen to shift downward for small values of x. 
Solar control of the D region extends to 70 km which, as we 
shall see in a future section, is indicative of large con- 
centrations of nitric oxide. The peak density of the D 
region resembles that of midlatitudes under quiet solar 
conditions. There is also an electron density minimum in 
the vicinity of the mesopause. This has been observed at 
midlatitudes and attributed to electron attachment to dust 
(Aikin et al., 1964). 

One interesting aspect of the comparison between pro 
files is the possible existence of a diurnal asymmetry, 
although this observation may be somewhat marginal within 
the accuracies estimated for the experiment. There is 
apparently a much larger change between x * 53° and 28° 
than between x * 28° and 48° below 100 km. The measurement 
for x ~ 48° occurred in the afternoon while the remaining 
two occurred in the morning. This asymmetry was also 
observed in medium frequency absorption measurements taken 
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during the day of the rocket flights at Colombo, Ceylon 
(Gnanalingam, private communication) . Similar ground- 
based observations have been reported previously (Rao et al., 
1962; Haug et al., 1970). 

SOLAR X-RAY PROPORTIONAL COUNTER EXPERIMENT 

Included in the payloads of rockets 14.425 and 14.372 

were X-ray proportional counters. The experiment consisted 

of an LND side window proportional counter tube whose gas 

filling was xenon plus 10% CH^ at 760 1 5 mm Hg. Each 

window diameter was 2.54 cm; the window material was 

-2 

beryllium of 2.54 x 10 cm thickness for 14.425 and 
5.08 x 10 ^cm for 14.372. 

The output of the proportional counter is fed into a 
charge sensitive amplifier, the output of which consists of 
pulses in the range 0 - 5 v dc. These pulses are stretched 
to a duration of 200 n sec with a rise and fall time of 
20 sec. Recovery time from the termination of the 
stretched pulse for a pulse of mean pulse height to 5% 
of the mean pulse height is 150m- sec. The resulting 
stretched pulses are telemetered to the ground over a 5 
KHz bandwidth telemetry channel. A pulse rate counter is 
included which telemeters a 0 - 5 volt signal that is 
proportional to the logarithm of the pulse repetition 
frequency. 

The unattenuated X-ray spectra obtained from the two 
payloads are shown in Figure 4. They compare favorably with 
Solrad 9 satellite ion chamber results below 5A. Kreplin 
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(private communication) has indicated the 1A to 8A integrated 
flux values for the local times of 14.425, 14.424 and 14.372 

o _ j n 

to be 1 x 10 , 1.3 x 10 and 9.0 x 10 ergs/cm sec, re- 
spectively. Kreplin (1961) has outlined the procedure for 
obtaining integrated flux from satellite X-ray ion chamber 
measurements such as those made on Solrad 9. Between 3 

0 Q 

and 8A, a 2 x 10 K gray body is assumed while below 3A, 

a 10 7 °K gray body is used. In each case an appropriate 

dillution factor is chosen. For comparison purposes two 

gray body curves are shown in Figure 4. Gray body I cor- 

6 o 

responds to a temperature of 2 x 10 K and a dillution 

factor of 5 x 10~ 15 . The integrated flux is 1.7 x 10” 3 ergs/em 2 

sec. Gray body II represents a temperature of 10 7 °K and a 

-21 

dillution factor of 1 x 10 . The integrated flux below 

—6 2 

3A is 2.5 x 10” ergs/cm sec, well below the threshold 
-5 2 

of 1 x 10 ergs/cm sec for the satellite 0.5 - 3A ion 
chamber. The proportional counter results are illustrated 
in Figure 4 and are seen to concur with the spectra estimated 
from Solrad 9 data. 

PRODUCTION FUNCTION 

In the absence of solar X-ray flares and energetic 
particle precipitation, the primary D region ionization 
source is solar Lyman alpha radiation acting on nitric 
oxide. Unfortunately the nitric oxide distribution has 
been one of the principal uncertainties in explaining the 
D region. Theoretical NO concentrations derived from 
photochemical equilibrium involving the reaction of ground 
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state nitrogen atoms with oxygen and nitric oxide are not 
large enough to account for the required ion-pair production 
function, Aikin et al. (1964). On the other hand NO dis- 
tributions derived from rocket measurements of resonance 
fluorescence of nitric oxide gamma bands (Pearce, 1969) 
would require drastic restructuring of the theory of electron 
loss processes and be inconsistent with observations of 
solar flare X-ray enhancements of D region ionization, 
(Soraayajulu and Aikin, 1968). However, recent observations 
of Meira (1971) have shown that in the D region the analysis 
of previous nitric oxide resonance fluorescence data did not 
consider the background due to Rayleigh scattering and so 
overestimated the nitric oxide concentration. 

Two similar rocket-borne spectrophotometer measurements 
were made by Meira at Wallops Island, Virginia in January, 
1969. The range of measured nitric oxide values is shown 
in Figure 5. An ion pair production function has been 
derived to correspond to the solar zenith angle for each 
rocket flight. It employs an NO mixing distribution passing 
through the lower limit of the error estimate of the best 
fit of Meira 's data. This value is further justified by 
the recent results of Rusch and Barth (1971) who have 
discovered a trend toward lower values of NO in the 
equatorial region. The resulting production function, 
q(NO) , for x * 28° is shown in Figure 6. 

An X-ray production function based on the rocket 
proportional counter spectra and satellite results is also 
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shown in Figure 0. Due to the hardness of the spectrum, 

X-rays are not a negligible contribution to the D region 
ionization rate and with cosmic rays constitute the major 
source for ionization of the major atmospheric constituents. 

The expression for the contribution due to cosmic rays is 

q - q Q n(sO/n Q 

where q Q is ion pair production function under standard 
conditions at the geomagnetic equator, i.e. 25 ion pairs/cm sec 
(Nicolet and Aikin, 1900). The number of particles in a 
standard atmosphere n Q , is 2,7 x 10 * cm . The resulting 
cosmic ray production function is also shown in Figure 6 
and merges with the X-ray ionization rate. 

It has been suggested by Hunten and McElroy (1968) that 

1 ° 

photoionization of 0 2 ( A g ) at X < 1118 A is an important 
source of D region ionization. Huffman et al. (1971) have 
pointed out that this would indeed be the case if there were 
no c ^2 in *k e mesosphere. Including the expected atmospheric 
COg absorption of the ionizing radiation leads to the ion 
pair production function shown in Figure 6. This is con- 
sistent with Huffmann's conclusions. The ionization of 
0g( 1 Ag) Is therefore, not a significant source of D region 
ionization. If anything it plays a minor role in the 
formation of the E region. 

The possibility that 0g can be raised to vibrationally 
excited states and subsequently ionized as suggested by 
Inn (1961, a, b) can also be discarded. It has been shown 
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by Hudson (1971) that such a sequence will lead to dis- 
sociation of the 0 2 molecule rather than ionization. 

The sources of ionization for the equatorial D region 
have been reviewed and estimated. We conclude that sources 
listed by Nicolet and Aikin (1960), i.e. ionization of 
nitric oxide by solar Lyman alpha, short wavelength solar 
X-rays and galactic cosmic rays, and a small contribution 

i 

due to 0 Q ( A ) ionization, are of primary importance for 

describing the D region ionosphere under quiet conditions. 

EFFECTIVE RECOMBINATION COEFFICIENT 

The term "effective recombination coefficient" is used 

2 

here to define the ratio, a m ? q./tt where q, is the 

i e i 

ion production rate for each ion species "i" and N 0 is the 
electron density. It is a useful parameter because it tells 
one if the chosen ion-pair production function is consistent 
with current ideas on the magnitude of the loss rate. 

The a*s which are derived from the ion pair production 
functions and the electron density profiles at x * 28° are 
shown in Figure 7 for two nitric oxide distributions at the 
lower and upper limits of Meira's NO measurement. Simi .ar 
results are obtained for x * 53° and 48°. They are quite 
significant when compared with recent measurements of 
Biondi (private communication) for the dissociative 

"f 

recombination coefficients of the water clusters ions 19 , 
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Measurements 


H 3 0*; 37*, H 3 0*.(H 2 0) 

; and 55 

- «3 

at 30C °K are modified 

by the 

term 

^lO* 

- 1 . 3 x 

10" 6 

*37* 

- 2.7 x 

CO 

1 

o 

r—4 

*55* 

- 4.6 X 

10" 6 

°30* 

“ 4 . 0 x 

10- 7 


l 2 '2 


TT J 


( bULu 
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The altitude variations for these coefficients are shown in 
Figure 7 based on a temperature model from Champion (1965) . 
Note that the effective recombination coefficients derived 
from our measurements at X « 28° bracket the values for a 
mean water cluster mass near 55* . It is shown in Goldberg 
and Aikin (1971) that 37*, 55* and heavier ions do indeed 
predominate below 80 km. This result presents a consistent 
picture between ionospheric and laboratory measuremfents 
for the first time since Aikin et al. (1964) pointed out 
the necessity for large recombination coefficients to 
accomodate the rocket measurements of nitric oxide. This 
indicates validity for the theory put forward by Reid (1970) 
and Landmark et al. (1970) that faster dissociative re- 
combination will occur if the species of ion is a water 
cluster ion than if it is 30*. In order to reconstruct the 
effective recombination coefficient it is necessary to take 
into account the proportions of the different species of ions. 
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CONCLUSION 


Three equatorial D region electron and positive ion 
density profiles have been presented for one day. Evalua- 
tion of the possible production sources leads to the conclusion 
that X-rays, photoionization of NO, and cosmic rays dominate 
ion production in the D region. X-ray fluxes are determined 
by rocket-borne proportional counter experiments aboard two 
of the rockets and are found to concur with the integrated 
total flux measured aboard the Solrad 9 satellite at that 
time . 

Comparison of the measured electron density profiles 
with the X-ray flux and recent measurements of nitric oxide 
provide effective recombination coefficients between 

A O _ C n 

4 x 1(T cm /sec and 1 x 10” cm /sec in accord with recent 
measurements of the electron-ion recombination coefficient 
for heavier hydrated cluster ions near mass 55 + . No ad- 
ditional mechanisms, e.g. large concentration of negative 
ions, are required to resolve the long-standing discrepancy 
between ion-pair production rates due to ionization of 
nitric oxide and laboratory measurements of loss processes. 

It was shown by Aikin (1966) that for recombination 

-6 3 -1 

coefficients as large as 10 cm sec , no diurnal asymmetry 
of electron density is^expected. Hence the apparent asym- 
metry about noon cannojtMae^explained at the present time. 

A large gradient in electron density above 85 km often 
observed at midlatitudes is also present at the equator. 
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TABLE LIST 


Table 1 


Summary of instrumentation, launch times and solar 
zenith angles for rockets 14.425, 14.424 and 14.372 
An X implies that the experiment was carried aboard 
the specified payload. 
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TABLE I 


Rocket Number 

14.425 

14.424 

14.372 

Solar Zenith Angle 

53.2° 

0 

00 

• 

48.3° 

Time of Launch (LMT) 

0827 

1017 

1509 

Electron density 
(radio propagation) 

X 

X 

X 

Electron density 
(nose tip probe) 



X 

Positive ion density 
(Gerdien condenser) 

X 

X 


Positive ion composition 
(quadrupole ion mass 
spectrometer) 

X 

X 


Solar, Lyman cr & 0„ density 
(ion chambers) 

X 

.. 

X 

X 

Solar X-rays 
(proportional counter) 

X 


X 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


Electron and positive ion density altitude 
distribution; la-14.425, lb-14.424, lc-14.372. 
Gerdien condenser for measuring positive ion 
number density. 

Comparison of electron density profiles for 
March 19, 1970 at Thumba , India. 

X-ray spectra obtained from 14.424 and 14.372. 
Gray body curves are described in the text. 
Comparison of nitric oxide distributions adapted 
from Meira (1971). Arrow on Barth’s curve refers 
to correction required for Rayleigh scattering. 

D region production function for a solar zonith 
angle of 28°. 

Effective recombination coefficient profiles 
derived from 14.424. Upper and lower limit 
curves refer to limits on Meira (1971) data. 

Shown for comparison are the dissociative 
recombination coefficients of HgO*, HgO^HgO, 
and H 3 0 + ’(H 2 0) 2 . 
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